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Abstract 

This paper introduces a systematic treatment of the linear theory of scalar grav- 
itational perturbations in the presence of a fully inhomogeneous magnetic field. The 
analysis is conducted both in the synchronous and in the conformally Newtonian gauges. 
The cosmological plasma is assumed to be composed of cold dark mattter, baryons, pho- 
tons, neutrinos. The problem of super-horizon initial conditions for the fluid variables 
of the various species and for the coupled system of Boltzmann-Einstein equations is 
discussed in the presence of an inhomogeneous magnetic field. The tight coupling ap- 
proximation for the Boltzmann hierarchy is extended to the case where gravitating 
magnetic fields are included. 



1 Introduction 



The impact of magnetic fields on the anisotropies of the Cosmic Microwave Background 
(CMB) has been a subject of active investigation since the early attempts of Zeldovich 
[1, 2], where it was proposed that all the anisotropy in the CMB could be generated by a 
magnetic field. Following the development of inflationary cosmology it is now clear that 
the explanation of the large-scale temperature anisotropies should be attributed to some 
adiabatic (or quasi-adiabatic) mode, which was present outside the horizon prior to the 
decoupling of radiation from matter. 

The theory of scalar gravitational fluctuations is a well developed subject and it is an 
essential tool setting the initial conditions for the evolution of CMB anisotropies. When 
initial conditions are set, usually deep within the radiation era (but after neutrino decou- 
pling taking place around 1 MeV), the dominant component of the plasma are photons and 
(effectively massless) neutrinos. The subdominant component of the plasma is formed by 
baryons, electrons, cold dark matter (CDM) particles. CDM particles are only coupled to 
gravitational interactions and behave like a perfect relativistic fluid. Also the photons, lep- 
tons and baryons, being tightly coupled through Thompson scattering, behave like perfect 
relativistic fluids. On the contrary, neutrinos are essentially collisionless and, therefore, do 
not really behave like a perfect fluid. Because of this physical difference, neutrinos should 
be described through the appropriate Boltzmann hierarchy of their phase-space distribu- 
tion. If the dark energy is paramertized in terms of a cosmological term, the gravitational 
fluctuations of this sector do not affect the problem of the initial conditions. 

If primordial fluctuations of the geometry are present outside the horizon before matter- 
radiation equality, they imprint also fluctuations in the density contrasts and peculiar ve- 
locities of the different species. The quantitative evolution of the various plasma quantities 
is determined by the fluctuations of the geometry by solving the coupled system of Ein- 
stein and fluid equations in the radiation-dominated epoch when the relevant modes are 
still outside the horizon. Following the classification scheme pioneered by Bardeen [3], the 
initial conditions may be classified into adiabatic or isocurvature modes. In the adiabatic 
case the total fluctuations in the entropy density of the CDM-photon-lepton-baryon fluid 
vanish at large distance scales which is compatible with a constant mode of curvature fluc- 
tuations outside the horizon. In the case of isocurvature modes, in the usual terminology, 
the fluctuations in the entropy density do not vanish; according to the linear analysis of 
the Einstein-fluid system, this is compatible, in some cases, with non-constant modes of 
the curvature fluctuations. In spite of this simplified classification, the situation may be 
more complicated since different isocurvature modes may be allowed for the different species 
present in the plasma. Moreover, it is also quite plausible to have a situation where the ini- 
tial conditions are, predominantly, of adiabatic nature, but with a subdominant isocurvature 
component. 

Although an excellent approximation early on, the tight-coupling assumption breaks 
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down at later times, when photons and baryons decouple. In spite of this caveat, a sys- 
tematic use of the tight-coupling expansion allows important analytical estimates of the 
produced CMB anisotropics [4]. 

In the theory of the CMB anisotropies, the problem of initial conditions is often presented 
in two complementary descriptions, namely the conformally Newtonian and the synchronous 
gauges. Both gauges are quite useful for different reasons (see also [6]). The conformally 
Newtonian gauge (often dubbed longitudinal) is effective in the discussion of the evolution 
of the fluctuations while they are still outside the horizon. However, in such a gauge, the 
discussion of the initial conditions for the CMB anisotropies may be rather complicated. For 
instance, there are physical isocurvature modes that are singular in the longitudinal gauge 
but not in the synchronous gauge. Furthermore, various codes needed for the numerical 
calculations of the CMB spectra are formulated within the synchronous coordinate system. 
There have been, in the past, controversies concerning the use of the synchronous gauge 
[5]. The main caveat was that gauge modes may appear in the synchronous gauge. These 
gauge modes, however, can be precisely disentangled from the physical ones [5], so that a 
pragmatic approach envisaged some time ago [6] has been to discuss the problem of initial 
conditions in both gauges in parallel. 

Suppose now that inhomogeneous magnetic fields are present. In this case two major 
effects can be envisaged: 

• inhomogeneous magnetic fields can gravitate, i.e. they can affect the perturbed Ein- 
stein equations, thus becoming sources of the fluctuations of the geometry; 

• inhomogeneous magnetic fields can impart different velocity gradients to the baryon- 
photon-lepton fluid, but not to the neutrinos and to the CDM components. 

Deep within the radiation epoch, large-scale magnetic fields can be treated, in first ap- 
proximation, as interacting with a single globally neutral fluid. In fact the two electrically 
charged species present in the problem are the electrons and the baryons which are in ther- 
mal equlibrium at a common temperature. The typical length scales of interest for the 
present discussion are the ones much larger than the typical magnetic diffusivity scale (set 
by the finite value of the conductivity) and also much larger than the Silk damping scale 
(set by the finite value of the thermal diffusivity scale). Furthermore, for this range of 
scales we are not interested in the propagation of high frequency electromagnetic waves 
whose specific analysis would clearly require a two- fluid plasma treatment. The one- fluid 
plasma description if often dubbed as magnetohydrodynamics (MHD). In MHD, the cur- 
rent density, the magnetic field and the electric field are all solenoidal and the Lorentz force 
affecting the baryon peculiar velocity modifies the standard tight-coupling approximation. 

By looking at the evolution of the brightness of the temperature fluctuations it is rather 
clear that, to zeroth order in the tight-coupling expansion, the monopole and dipole equa- 
tions are modified. It is known [4] that to zeroth order in the tight coupling expansion 
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the CMB is not linearly polarized since, to this order, photons and baryons are so strongly 
coupled that the photon distribution is isotropic in the baryon rest frame. The photon dis- 
tribution being isotropic, Thompson scattering does not polarize the CMB. To first order in 
the tight-coupling expansion the polarization of the CMB is proportional to the quadrupole 
of the photon distribution. The quadrupole in the temperature fluctuations is, in turn, 
produced by the free streaming of the dipole between collisions. If a Lorentz force term is 
present, the evolution of the monopole acquires a further source term and, as a consequence, 
all the higher orders in the tight coupling expansion are modified. 

Various analyses of the possible impact of magnetic fields on standard CMB physics 1 
have been proposed up to now (see, for instance, [7] for a short review on this specific 
subject). In [8] the analysis of scalar fluctuations of the geometry is made under the 
assumption that the Lorentz force vanishes. In MHD this assumption translates into a 
specific class of magnetic field configurations. However, even in the case of a force-free 
configuration the specific initial conditions used by the authors in order to compute the 
CMB anisotropies are not specifically discussed, i.e. no solution for super-horizon sized 
fluctuations including large-scale magnetic fields, is presented. Furthermore, no specific 
investigation of the impact of gravitating magnetic fields on the tight-coupling expansion 
has been attempted. In [9] an interesting analysis of the impact of large-scale magnetic 
fields on CMB anisotropies has been performed, mainly for the vector and tensor modes of 
the geometry. The scalar fluctuations have not been discussed. In [10] attention is always 
paid to vector and tensor modes and the problem of initial conditions is more accurately 
specified. It is clear that without a precise specification of the initial conditions for the 
CMB anisotropies the numerical analysis is rather difficult to implement since it is unclear 
what should be, for instance, the initial values of the lowest multipoles of the phase-space 
distribution for the various species. The need of an accurate analysis of the problem of initial 
conditions can be also understood by observing that the system of perturbed equations, even 
if linear, has many unknowns so different solutions, describing different physical situation 
may be possible. 

The effect of large-scale magnetic fields on CMB anisotropies may also be discussed in 
the case when the magnetic field is not completely inhomogeneous. There could be situations 
in which there is an homogeneous component of the magnetic field (for instance along the 
z axis. This case is similar to the one originally investigated in [1, 2]. In this framework 
bounds on such a magnetic field can be derived [11]. It could also happen that the uniform 
magnetic field supports inhomogeneities in the bulk velocity field affecting, consequently, 
the CMB anisotropies. One example in this direction are the Alfven waves, whose effects 
and implications have been analysed under different approximations by various authors 
[12, 13, 14, 15]. The simplification of having a uniform magnetic field also allowed the 
analysis of possible Faraday rotation effect of the CMB polarization plane [16, 17, 18, 19]. 

x By standard CMB physics we mean that a constant mode of adiabatic (or quasi-adiabatic) fluctuations 
is assumed to be present outside the horizon prior to equality. 
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The motivation for including a uniform magnetic field relies on the simplicity of the 
configuration. However, the magnetic fields produced in the early Universe, as far as we 
understand the problem, are unlikely to be produced in a perfectly uniform configuration. 
They are most likely to be produced either from the amplification of vacuum fluctuations 
of some primordial gauge field or from some phase transitions. In all these cases magnetic 
fields are expected to be fully inhomogeneous (see, for instance, [21] for a recent review on 
the role of large-scale magnetic fields in cosmology and astrophysics) . 

The purpose of the present study is to give a systematic analysis of the theory of scalar 
gravitational fluctuations in the presence of a fully inhomogeneous magnetic field and in the 
context of the conventional Bardeen formalism. This analysis is mandatory both to derive 
accurate bounds on large-scale magnetic fields from CMB physics and to bring the study 
of magnetized initial conditions to the same standard as in the non-magnetized case. To 
the best of our knowledge this problem did not receive specific attention in the literature 
concerning the linearized theory of scalar fluctuations in the presence of magnetic fields. The 
results reported in the present paper are general in the sense that no specific configuration 
of the magnetic field is assumed. The only assumption, as previously stressed, is that the 
magnetic field is fully inhomogeneous. 

It is important to mention that various studies discussing the evolution of large-scale 
magnetic field exist in the covariant approach which is somehow complementary to the one 
adopted in the present paper. For a full discussion of the problem see [22, 23] and references 
therein. 

The plan of the present paper is the following. In Section 2 the evolution equations 
for gravitating inhomogeneous magnetic fields will be derived in the MHD approximation 
and in the context of the conformally Newtonian gauge. In Section 3 the initial conditions 
for the CMB anisotropies will be analysed. Solutions for super-horizon-sized fluctuations 
will be presented in different cases. Section 4 is devoted to the analysis of the synchronous 
gauge description. In Section 5 the tight-coupling expansion will be revisited in the presence 
of inhomogeneous magnetic fields. Section 6 contains some concluding remarks. Various 
useful technical results have been collected in the appendix even though, in some cases, they 
were already present in the literature, but within a different set of conventions or within a 
different context. 

2 Gravitating inhomogeneous magnetic fields 

The treatment of magnetic fields in a globally neutral plasma differs considerably according 
to the specific range of scales relevant for the problem under discussion. For instance, the 
electromagnetic branch of the spectrum of plasma excitations, in a cold or warm magnetized 
plasma, can only be partially addressed within a two-fluid description [24, 25]. Partially 
means, in the present context, that the dispersion relations can be obtained but the typical 
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damping scales require the solution of the full kinetic system of equations for the different 
species (the so-called Vlasov-Landau approach). 

The phenomena occurring over much smaller frequencies and over large length scales can 
be described on the basis of our knowledge of terrestrial plasma, by an appropriate one-fluid 
description, which is provided by MHD [26]. The main idea behind the MHD approximation 
is, in short, the following. Starting from a two-fluid description (for instance electron and 
baryon fluids) it is possible to define appropriate one-fluid variables. For instance the bulk 
velocity of the plasma appearing in the MHD description is, in flat space-time, 

_ m h v h + m e v e 

v = . (2.1) 

m h + m e 

which is, in practice, the centre-of-mass velocity of the baryon-electron system. In the spe- 
cific case touched by the present discussion, both baryons and electrons are non-relativistic 
and in thermal equilibrium T e b- 

The reduced MHD description [26] has been employed in the analysis of different prob- 
lems arising in connection with large-scale magnetic fields. For instance, the evolution of 
these fields in curved space-times is normally discussed within a MHD approach [27, 28, 29] 
(see also for a review [30, 31]). More formal discussions on MHD and on two- fluid de- 
scriptions in curved space-times can be found in [32, 33, 34]. Finally, the reduced MHD 
description has been used in the analysis of the implications of large-scale magnetic fields 
on structure formation [35, 36]. 



2.1 The perturbed system of Einstein equations 

Consider now the Friedmann-Robertson- Walker (FRW) line element written in the confor- 
mally flat case 

ds 2 = a 2 (i])[dr] 2 -dx 2 ], (2.2) 

where r\ is the conformal time coordinate. Since the magnetic fields are fully inhomogeneous, 
they will not affect the homogeneous background whose equations are, in the spatially flat 
case: 

H^X^E^ (2-3) 

6 A 

U 2 -H' = 4TrGa 2 pY,tt x {l + w x ), (2.4) 

A 

p'x + 3H(Px+Px) = 0, (2.5) 

where H = (In a)' and the prime denotes a derivation with respect to rj. In Eqs. (2.3)-(2.5) 
the summation index A refers to each component of the plasma, i.e. baryons, photons, 
neutrinos and CDM particles; w\ = p\/ p\ and Q,\ are, respectively, the barotropic index 
for the given species A and the energy density (in critical units) for each component of the 
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fluid. In view of the subsequent discussions, it is appropriate to define here also the usual 
Hubble parameter, i.e. H = a /a, where the overdot denotes a derivation with respect to 
the cosmic time coordinate t, which is related to r\ by the differential relation dr\a(j]) = dt. 
In terms of these definitions we also have H = aH. 

Let us now consider the fluctuations of the homogeneous FRW metric (2.2). In this 
section the conformally Newtonian coordinate system will be discussed and, therefore, the 
metric (2.2) can be consistently perturbed in terms of the two longitudinal scalar degrees 
of freedom (j) and ip 

$900 = 2a 2 '■(/), bgij = 2a 2 ip5 ij . (2.6) 

Defining the fluctuations in the Ricci tensor and in the Ricci scalar as 5R^ and 5R, the 
perturbed Einstein equations, 

SR; - ±5»5R = 8nGa 2 (5T» + 6t»), (2.7) 

relate the fluctuations of the geometry to the fluctuations of the matter sources, ST^, and 
to the fluctuations of the electromagnetic energy-momentum tensor, 5t". 

In explicit terms the (00), (Oi) and (ij) components of Eq. (2.7) lead, respectively, to 

V 2 Y> - m{H(j) + ip') = 47rGa 2 [<5T ° + 6t§], (2.8) 
-d\U(l) + ^) = 4mGa 2 {5T l Q + 5t 1 ), (2.9) 

si 



where 



4>" + HW + 4>') + (2H' + H 2 )<p + ^V 2 (0 - V) 
-\did j {<t> - i/>) = -^Ga 2 [5Ti + <5r/], (2.10) 



^ = ^(E 2 + B 2 ), (2.11) 
Sri = x B, (2.12) 



<5r/ 1 



EiE* + BiB* - \[B 2 + E 2 )S{ 



(2.13) 



4vra 4 

<5T ° = ]Tp A <5 A , (2.14) 

A 

5Ti = -'Ew x p x 6 x 6> i +V i , (2.15) 

A 

STt = y Z(l + w x )p x v\, (2.16) 

A 

where 5 X = 5p x /p x and v\ are, respectively, the density contrast and the peculiar velocity 
for each particle species; Y? i is the traceless component of the energy-momentum tensor of 
the fluid sources, i.e. 

A = m~ (2.17) 
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According to the usual notation the anisotropic stress can also be written, in Fourier space, 
as 

Q = d 3 d% = -k 2 ]T a x (l + w x )px, (2.18) 

A 

where a\ denotes the fraction contributed by each species A to the total anisotropic stress. 
Notice, as it will be more extensively discussed later, that the dominant source of anisotropic 
stress comes, for temperatures below the MeV, from the neutrinos. 

Finally in Eqs. (2.11)-(2.13) the components of the energy-momentum tensor have been 
written directly using the rescaled electric and magnetic fields E and B, which can also be 
expressed in terms of the corresponding flat-space fields as 

E(r],x) = a 2 (r])i(r],x), B(rj,x) = a 2 (r])B(r], x). (2.19) 

In the following the evolution of the different quantities appearing in the right-hand side of 
Einstein equations will be discussed in detail. 



2.2 Electromagnetic fields and baryons 

The Maxwell equations can be studied using the rescaled fields proposed in Eq. (2.19) 

B' = — V x E, (2.20) 

V x B = Air J + E', (2.21) 

V ■ E = 4irp q , (2.22) 

V ■ B = 0, (2.23) 

where p q is the charge density (p q ~ in a globally neutral plasma); J = a 3 j is the Ohmic 
current density 

J = a(E + vxB). (2.24) 

In Eq. (2.24), a is the conductivity that is related to the flat-space conductivity as a = aa c . 
Notice that in the reduced MHD description the Ohm law follows, in the reduced MHD 
approach, by taking the difference of the equations describing the momentum conservation 
for electrons and for ions. This procedure is rather tricky and, in principle other terms 
may appear in the generalized Ohm law such as the thermoelectric term and the Hall term 
[24, 21]. These will be neglected here. Finally, as pointed out after Eq. (2.1), the bulk 
velocity field appearing in Eq. (2.24) can be identified with the velocity of the baryons. 

The evolution for the baryon velocity field can be derived by perturbing the covariant 
conservation equation (including the appropriate electromagnetic contribution), i.e. 

d,5TC + 6T» a T£ a + +T; a ST£ a + ST^ + T^ST^ - F™j a = 0, (2.25) 

where T^" and 8Tt v are, respectively, the energy-momentum of the baryons and its first- 
order fluctuation; F va is the Maxwell field strength and j a is the current density. An 
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analogous notation is used for the Christoffel symbols computed on the background (denoted 
by an overline) and for their first-order fluctuations. Thus, defining 6>b = V--Ub, the evolution 
equation for the peculiar velocity of the baryons and for the density contrast reads: 

e' h = -ne h - c 2 s v 2 5 h - v 2 </> + ^ ane x e a T (e, - e b ) + v ^ JxB \ (2 . 2 6) 

5' h = W - 9 h . (2.27) 

where x e is the ionization fraction of the plasma, <tt the Thompson cross section; Q 7 Q b 
are, respectively, the energy densities, in critical units, for photons and baryons. At this 
point we can also define the differential optical depth for Thompson scattering: 

t 1 = x e n e aTa(r)). (2.28) 



Equations (2.20)-(2.23), together with Eqs. (2.24) and (2.26), can be studied in the 
MHD approximation. Since the plasma is globally neutral (p q ~ 0), the charge density of 
the electrons will be exactly compensated by the baryons and, according to Eq. (2.22), the 
electric field will be solenoidal. Moreover, for sufficiently small frequencies, the displacement 
current can be neglected in Eq. (2.21). This implies that also the Ohmic current will be 
solenoidal, i.e. 

V-J = 0, J=^-VxB. (2.29) 

As a result, the Ohmic electric field vanishes exactly in the limit of vanishing conductivity, 
i.e. since 

- J V x B , 

E = , (2.30) 

a a 

the electric fields appearing in Eqs. (2.11) and (2.13) lead to a contribution that is sup- 
pressed as 1/a 2 . Similarly, the electric field appearing in Eq. (2.12) leads to a contribution 
that is suppressed as 1/a. 

Recalling, as previously defined, that a = a c a, for temperatures T < 1 MeV the con- 
ductivity can be expressed as [24, 21] 



J_(T±) 1/2 

«cm \m c J 



(2.31) 



where T e b, as previously defined, is the common temperature of electrons and baryons. The 
finite value of the conductivity sets a typical diffusion scale for the evolution of the magnetic 
fields which is typically of the order of 



72.2 \MpJ 



La(T) = ^-^-l^-\ ^ L H (T), (2.32) 



where L a is the physical diffusion scale, Lh — H 1 is the Hubble radius and the number 
of relativistic degrees of freedom. Since, in our problem eV < T eb < MeV, L a <C Lh which 
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means that the finite value of the conductivity only affects scales which are much smaller 
than the Hubble radius. 



Similar, but quantitatively different, remarks can be made for the thermal diffusivity 
[21]. For temperature smaller than the temperature of neutrino decoupling, photons are the 
most efficient source of momentum transfer and the thermal diffusivity scale can be written 

as 

. 1.03 x ^f^VV^y'V^V' 4 ,/^, (2.33, 



Lh(T) ' ' V0.02 7 V0.5/ V 10.75/ V T * 

where, according to previous definitions (see Eq. (2.26), x e is the ionization fraction, 
the energy density of baryons in critical units. 

2.3 Cold dark matter 

The evolution of the CDM particles follows from the fluctuations of the covariant conserva- 
tion equation whose first-order fluctuation leads, in Fourier space, to 

e' c + ne c = k 2 ^, (2.34) 

s' c = 3V>' - e c . (2.35) 
where, following the previous notation 6 C = diV l c and <5 C = 5p c /p c . 



2.4 Photons and massless neutrinos 

For the neutrinos the evolution equations follow only partially from the energy-momentum 
conservation. In fact neutrinos at early times after neutrino decoupling (occurring around 
1 MeV) obey the collisionless Boltzmann equation (see the appendix for a more detailed 
discussion). To describe neutrinos during and after horizon crossing requires a Boltzmann 
hierarchy for 8 V: d v and for the higher multipole moments, i.e. t > 2, of the neutrino phase- 
space density T v i. With this caveat in mind, after neutrino decoupling, at temperatures of 
about 1 MeV, massless neutrinos obey, in Fourier space, the following set of equations (see 
the appendix for a swift derivation within the metric conventions adopted in Eqs. (2.2) and 
(2.6)) 

S' v = ~9 v + 4:iJ/, (2.36) 
k 2 

e' u = —5 U - k 2 a u + k 2 (j), (2.37) 

oi = ^=e v - ^fcTvs, (2.38) 
15 10 

where o v = ^2/2 is the quadrupole moment of the (perturbed) neutrino phase-space 
distribution and, as introduced above, T v i is the £-ih multipole. 
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The photon and neutrino evolution equations differ in the presence of an anisotropic 
stress term o v . Because of their frequent scattering by charged leptons and baryons, photons 
are unable, at early times, to develop a quadrupole moment in their velocity distribution. 
As a consequence, when a mode enters the horizon, the photons behave as a perfect fluid, 
while the neutrinos free stream, creating inhomogeneities in the energy density, pressure and 
momentum density. This is the physical reason why the anisotropic stress term appearing 
in Eq. (2.15) is effectively dominated, at early times, by the neutrino contribution. Also 
for the photons it is more appropriate to discuss the full Boltzmann hierarchy, especially in 
the light of the subsequent analysis of the tight-coupling approximation. This analysis will 
be summarized in the appendix. In the fluid approximation the evolution equations of the 
photons are given by the covariant conservation of the energy-momentum tensor, i.e. 

<S; = ~0 7 + V, (2.39) 
k 2 

= — <5 7 + k 2 4> + an c a T (e h - 7 ). (2.40) 

At early times the characteristic time for the synchronization of the photon and baryon 
velocities is Tb 7 ~ (n c cjT) _1 , which is small with respect to the expansion time and to the 
oscillation period. Thus the tight-coupling between photons and baryons means that, in the 
first approximation, <tt — ► °o and # 7 ~ #b- A more detailed analysis of the tight coupling 
expansion is postponed to Section 5. 



3 Large-scale solutions 

The full system of scalar fluctuations in the longitudinal gauge will now be consistently 
discussed and solved. Prom Eq. (2.8) we obtain, in Fourier space 2 and using Eq. (2.3): 

-m{H4> + y/) - fcV = ^n 2 [(R u 5 u + (i - i?„)<5 7 ) + n B (k) + n h s h + n c s c ], (3.1) 

where, for N u species of massless neutrinos, 

7 / 4 \ 4 /3 ft 

R= 8 Nu \n) ' Ru = TTr' R ^ = l ~ R ^ ( 3 - 2 ) 

so that R v and i? 7 represent the fractional contributions of photons and neutrinos to the 
total density at early times deep within the radiation-dominated epoch. 

In Eq. (3.1) the contribution of the magnetic energy density has been parametrized, in 
Fourier space, as 

fl B (k, V ) = ^ = — Lj / Bi{\p- k\)B\p). (3.3) 

2 In order to avoid possible confusions with subscripts referring to the different species of the plasma we 
will avoid to introducing a further subscript labelling the Fourier mode. 
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The appearance of convolutions in the expressions of the magnetic energy density is a 
direct consequence of the absence of a uniform component of the magnetic field whose 
background contribution, as repeatedly stressed, is vanishing. Notice also that p appearing 
in Eq. (3.3) is the energy density of the dominant component of the fluid, so, for instance, 
in the radiation-dominated epoch p = p r , a A p r ~ constant and Qb — constant. 

From Eq. (2.9), using the same notation as in Eq. (3.1), the following simplified equation 
can be obtained for the momentum constraint 



k 2 (H<f> + ^') = hi 2 



-(R^O^ + (1 — Ru)9-y) + + ^b^b + & C Q C 



(3.4) 



where FQ(k,r)) is the Fourier transform of the generalized Lorentz force, i.e. 

V • [(V x B) x B] = J d 3 kF B (k)e in - s , (3.5) 
where F&(k) is given by the following convolution 

d 3 p[(k.p)B J (p)B^\k-^)-(k l - Pi )B J (\k-^)B J (\k-p\)B^p)]. (3.6) 



Notice that a consistent treatment of fully inhomogeneous magnetic fields implies that 

08 <<: W <<: 11 (3 ' 7) 

in order not to over-close the Universe. 

In the ideal MHD limit 3 , i.e. a — > oo, the contribution of Fs(k) disappears from Eq. 
(3.4). This occurrence is not general: it will be shown in a moment that, even in the ideal 
MHD limit, the generalized Lorentz force does contribute both to the baryon evolution 
equation and to the anisotropic stress. In fact the evolution equation for the baryons 
becomes, in Fourier space and with the conventions used so far 

e> h + He h = k 2 <t> + ^^. (3.8) 

By taking the trace of Eq. (2.10) it is easy to obtain 

< + (2^ + <j>')H + (2H' + H 2 )<j) - y ((j> - VO = \n 2 [{R v 5 v + (1 - i?,))<5 7 + fi B (fc)]> (3-9) 

where the contribution of the other species vanishes since both baryons and CDM particles 
have a vanishing barotropic index, i.e. = w c = 0. 

3 Notice that, so far, different quantities have been denoted by a (with different indices), namely the 
Thompson cross section (i.e. or), the conductivity (i.e. a) and the fractional contribution of each species A 
to the anisotropic stress (i.e. ax). Since all the quantities are properly defined and used, there should be no 
confusion. We felt that it would be even more confusing if we changed the usually accepted notations for 
these three important quantities. 
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By taking now the difference between Eqs. (2.10) and (3.9), we can get the equation for 
the anisotropic stress, i.e. for the perturbed components i ^ j of the Einstein equations: 

di&W - V) - lv 2 (4> - V) = 8nGa 2 p^l + -^— A [b^ - ^ 2 ^'] J. (3.10) 



Applying now the differential operator <9j<9 1 to Eq. (3.10), the following relation can be 
obtained 

i i 

(3.11) 



V 4 (0- V) = 12^Ga 2 
where, following the conventions of Eq. (2.18) 



a = = -k 2 a„(p„ + ft,) = - |*V (3.12) 

In order to write Eq. (3.12) in an explicit form it has been assumed that the only anisotropic 
stress arising from the fluid sources is the one related to the quadrupole moment of the 
neutrino phase-space density. Going to Fourier space, Eq. (3.11) becomes 



4 

-g0V - ctb(&) 



(3.13) 



having defined 



_n B (k) F B (k) 

The quantity cjb can be interpreted as the anisotropic stress arising thanks to the inhomo- 
geneous magnetic field and, in the force-free case, <tb = ^b/3. 

In order to discuss the initial conditions for the CMB anisotropics, the usual procedure 
is to solve the system of equations for the fluctuations while the Fourier modes of the 
various fields are outside the horizon, i.e. krj < 1. Different sorts of initial conditions are, 
in principle, possible. An interesting situation arises when both ip and (f> are, in the first 
approximation, constant outside the horizon during the radiation-dominated epoch. In this 
case, and in the absence of magnetic fields, it is known that the system of scalar metric 
perturbations admits an adiabatic mode. Adiabatic means, in the present context, that the 
total fluctuations in the entropy density of the CDM-baryon-radiation fluid vanish at large 
length scales. This condition implies specific relations between the various density contrasts 
deep within the radiation-dominated epoch, i.e. 



S c ^ 6 h ~ - A 5 V ~ jj<V (3.15) 



In the following, after reviewing the specific analytic form of the adiabatic solution, new 
solutions, arising in the presence of an inhomogeneous magnetic field, will be illustrated. 
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3.1 The standard adiabatic scenario 



If the magnetic fields are absent, i.e. F& = Qb = the standard system of equations is 
recovered and the initial conditions for the various fields can be found by solving, simul- 
taneously, the baryon-photon system (i.e. Eqs. (2.26)— (2.27) and Eqs. (2.39)— (2.40)), the 
CDM and neutrinos equations (i.e. Eqs. (2.34)-(2.35) and Eqs. (2.36)-(2.38)) together 
with the perturbed Einstein equations whose specific form has been derived above from Eq. 
(3.1) to Eq. (3.13). In the following, only the final result will be given: 

t t 3 (525 + 188*, + 16*^ 2 2 

6h = 6c = "2*° 60 (25 + 2*,) V ' (3 - 16) 

= -20o - (^ + 188^ + 16^) 

7 ^ u 45 (25 + 2*,) Y 1 ' v ; 

- (75 +UR u -8 Rl) , 2 2 

* = ^ ~ 90 (25 + 2*,) *° k V ' (3 - 18) 

—r- ( 2 \ , (75 + 79*, + 8*;*) , l2 2 

*= (l+ g*)*. - 90 (25 + 2^) (3.19) 

-j 00,2 (65 + 16*,) ,, 43 , . 

^ = Y fc?? - 36 (25 + 2*,/ °^' (3 - 20) 

* b = - (^+14^- 8^) 

2 ' 360 (25 + 2 *,) ru 1 ' v ; 

- _ 0o ,o (75+14 *,-8*3) 43 

^ c " T M " 360 (25 + 2*,) 77 ' (3 - 22) 



_ 00,2 2 (65 + 16*,) 44 . . 

= 15* " - 540 (25 + 2*,) ^ " ' (3 - 24) 

where the overline in the various quantities has been introduced for future notational con- 
venience. Concerning this solution, a few comments are in order: 

• while the leading order of the solution is usually quoted (see for instance [6, 37]), the 
second is useful in the present case, in order to check what happens at the horizon 
crossing of a given mode (i.e. krj ~ 1) and in order to compare with the situation 
where the magnetic field is present; 

• while to leading order all the peculiar velocities of the plasma coincide, to next-to- 
leading order they are different; 

• the adiabaticity condition given in Eq. (3.15) holds not only to leading order but also 
to next-to- leading order in \krj\ 2 . 



(25 + 8*,) x , 2 2 
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3.2 Force-free case 



Consider now the case when the magnetic field is force-free, and suppose, as previously 
discussed, that the solution at large scale is, in the first approximation, a constant mode 
for the longitudinal fluctuations of the metric. Then, the solution can be parametrized as 

0~0o + ^fcy, ^ ~ Vo + CV^V ■ (3.25) 

We are interested in the solution of the system deep in the radiation epoch where a(r]) ~ rj 
and the dominant component of the fluid energy-momentum tensor consist of photons and 
neutrinos, i.e. Ob <C 1 and Q c <C 1. From Eq. (3.1), the density contrasts for photons and 
neutrinos can be written as 

5 7 = -200 - n B + A^r] 2 , (3.26) 
K = -20o - n B + A u k 2 V 2 , (3.27) 

where the constants Ay and A v are related to and by the following algebraic equation 

+ 2CV + y = - X -\R V A V + (1 - R U )A 7 }. (3.28) 

Inserting Eqs. (3.25) together with Eq. (3.26) into Eqs. (2.39) and (2.40), the following 
parametrization for the peculiar velocity of the photons can be derived 



k 2 



2 ~ 4 



■q + L> 7 /cV (3.29) 



together with two algebraic relations determining the constants 

3D 7 = ^ + CV, (3.30) 

A, = -^ + ^ + iC,. (3.31) 

With a similar procedure the baryon and CDM fluids can be analysed: from Eqs. (2.26)- 
(2.27) and Eqs. (2.34)-(2.35), the solution can be written as follows 

0b = ^M + A>fcy, (3.32) 



and 



^p-^ + D^W, (3-33) 



S h = -^(20 o + n B ) + A h k 2 r l 2 , (3.34) 
6 c = -^(20 o + n B )+A c k 2 r] 2 . (3.35) 
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The constants appearing in Eqs. (3.32)-(3.35) should satisfy 

A - qr ^° 
4Aj,c = Cs, 



(3.36) 
(3.37) 



where the notation A^ c and c means that the relations hold, separately, for the baryon 
and CDM integration constants. 

The anisotropic stress a u appears both in Eqs. (2.37) and (2.38) and in Eq. (3.13), 
determining the difference between the two longitudinal fluctuations of the metric. Hence, 
from Eqs. (3.13) and (2.36)-(2.37) the neutrinos fluid variables are determined to be 



k 2 7] 2 D v 4 4 

<Tv = —=—{1p0 - 4>0) - — ^- + — K 7] , 



k 



4R U 15 



f^B R u — 1 



4^,3 



subject to the constraints 



2 4 R u 
5 V = —(2<?!>o + ^b) = A u k 2 rj 2 , 



o n A v _ n 00 ^B Ru ~ 1 
3D U - — _C -- + W -^— . 



— C,/ 



The lowest order of Eq. (2.38) also implies 



V'o 



7 A, 



b \l + -R v 



{Ru - l). 



(3.38) 

(3.39) 
(3.40) 

(3.41) 
(3.42) 
(3.43) 

(3.44) 



Up to now all the evolution equations for the fluid variables have been consistently solved 
in terms of a set of unknown constants satisfying a set of (algebraic) consistency conditions. 
From the remaining two Einstein equations (3.4) and (3.9), two further constraints on the 
various constants can be obtained, namely 



CV + 2CV, = 2D U R V + 2(1 - R U )D 7 , 
6C^ + C<f> - 



2 £1 1 

—Ru4>o + -^(Ru-l) + ~[A V R V + (1 - R u )A y ]. 



(3.45) 
(3.46) 



Since two of the algebraic equations are not independent, the system of ten algebraic equa- 
tions in ten unknowns becomes 



A h = A c = - 



69 - 61 R v 



60 (25 + 2^) 



Rl _ 525 + 188 R v + 16 R 2 , 



60 (25 + 2i2„) 
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A v = - 

Ccj, = — 

C^jj = — 
Dj = — 
D u = - 
D h = D c 



237 + 152^ + 16^ 
90 (25 + 2i2„) 



2 (525 + 188^ + 16^) 



45 (25 + 2ik) 



375 - 207 R u - 152 R 2 V - 16 R? v (525 + 188 R v + 16 R 2 ) 



90 R v (25 + 2R V 



Q B + 



45 (25 + 2i2„) 



6-8R u + 2R 2 75 + 14^-8^ 

L 45 (25 + 2 " ■ ■ - — 

69 - 61 Ik - 



90 (25 + 2 



180 (25 + 2i? 
(25 + 8 Ik 



8i?2 (75 + 79^ + 8i?2) ' 

"B H , z „ N 00 



</>0 



90 (25 + 2ik) 
(7 + 81k) 



.20 (25 + 2 Ik) 40 (25 + 21? 
45 - 29 R v - 16 Rl n (65 + 16 R, 



72 R v (25 + 2i?^ 



36 25 + 2 Ik 



Q-8R U + 2R 2 U ^ 75 + 14^-8^ 
180 (25 + 2 Ik) B + 360 (25 + 2 Ik) 



(3.47) 



Recalling the form of the standard solution reported in Eqs. (3.16)— (3.24), and using 
the results obtained so far, the full solution can be written as 



69 - 61 R - 8 Rl 



5» 



n B + 

8„ — f^B - 



60 (25 + 2 Ik 
237 + 152 Ik + 16 Rl~] 



Q B k rj 



90 (25 + 2 Ik) 
375 - 207 R v - 152 R 



n B k 2 v 2 , 



16 US 



901k (25 + 2 Ik) 
6 - 8 Ik + 2 i?2 



2^2 



Q B k Z 7] 



1p = ijj — 



45 (25 + 2 Ik) ^ B 
69-61^-8^2 



; 2 2 



180 (25 + 2 Ik) 



; 2 2 



2 (7 + 81k) 
4 ^ + 40 (25 + 212) 



5^B 



45 - 29^- 16 i?2 
72 Ik (25 + 2 ik) 



9 -a _^^Lk 2 n 



n B 



4^3 



#b — &c — 0b — 



6-8R„ + 2R 2 y 



180 (25 + 2 Ik) 



(3.48) 
(3.49) 
(3.50) 
(3.51) 
(3.52) 
(3.53) 
(3.54) 
(3.55) 



where the overline in the various quantities in the left-hand side denote the adiabatic solution 
derived above. 

Concerning these results, the following comments are in order: 



• unlike the standard case, the difference between the two constant modes of the longi- 
tudinal fluctuations of the metric is also determined by the magnetic energy density; 
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in the limit fie — ► 0, the standard expression (obtained in the previous subsection) is 
recovered; 



• the adiabaticity condition is enforced to lowest order in krj, but it is violated at the 
next-to-leading order; 

• the force-free approximation neglects, by assumption, the Lorentz force and this is the 
rationale for the equality, for instance, between 6^ and 9 C : in the absence of Lorentz 
force the magnetic field does not impart gradients to the baryonic component of the 
fluid; 

• the presence of magnetic fields introduces a further source of anisotropic stress as it 
can be appreciated from the chain of (first-order) relations ^ ip— 4>— Qb(Ri/ — l)/5 
which follows from Eq. (3.44). 

Since the adiabaticity condition is only realized to lowest order in kr) it is legitimate to name 
the solution presented in Eqs. (3.48)— (3.55) quasi- adiabatic. The corrections to adiabaticity 
are of order fl^k 2 !] 2 , so they are suppressed outside the horizon and also, by virtue of Eq. 
(3.7), by J1b- Notice that three regimes emerge naturally. The quasi- adiabatic regime where 
^b(&) < <j>o{k), the isocurvature regime Ob(&) > 4>o{k) and the fully adiabatic regime, i.e. 



The force-free limit may be (approximately) realized in nature in some specific (but 
rather interesting) cases. For instance, it could happen that magnetic fields are generated 
in a so-called maximally helical configuration. Recall, in fact, that given a magnetic field 
configuration, in a MHD description, it is always possible to define the related magnetic 
helicity [26], which is 



Now, at finite conductivity, the evolution of the magnetic helicity can be written as 



where the expression appearing in the right-hand side is called magnetic gyrotropy [26] and 
measures the number of contact points in the magnetic field lines. Clearly, if the magnetic 
gyrotropy is maximal, the orthogonal combination, namely the Lorentz force V x B x B , 
will be minimal or even zero. Maximally helical configurations can be very relevant in the 
development of MHD turbulence [26] and have been extensively studied in a cosmological 
context, (see also [21] for a general discussion and [38] for a recent paper). In particular it 
should be noticed that helical configurations can be produced at the electroweak time by 
coupling the hypercharge field to a pseudo-scalar [39] (see also [40]). 



MA:) < Mk). 




(3.56) 




(3.57) 
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3.3 Contribution of the Lorentz force 



Consider now the case when the magnetic field is not force-free. The major difference with 
the force-free case is that the evolution equation for the peculiar velocity of the baryons 
receives a contribution from the magnetic field gradient. Because of the tightly coupled 
dynamics of baryons and photons, also the photon peculiar velocity will be modified. Since 
photons are related to neutrinos by the Einstein equations, the whole solution will be 
modified. 

Equation (2.26) can be written, after multiplying both sides by the scale factor a(rj), as 

(a9 h )> = k 2 <j) a + — ^— , (3.58) 

where it has been assumed that 4> has, to lowest order, a constant solution denoted, as 
usual, by 4>G- The second term in the right-hand side of Eq. (3.58) is then constant and 
direct integration gives 

From Eq. (2.27), using Eq. (3.59), the density contrast can be determined up to a constant, 
which is fixed from the adiabaticity condition: 

5 h = - 3 -4>o ~ |0b - (3-60) 

The Hamiltonian constraint of Eq. (3.1) and the equation for the peculiar velocity of the 
photons imply 

8 V ~ S 1 ~ -20o - fi B , (3-61) 

The CDM equations (2.34) and (2.35) will then turn out to be 

<5 C ~-20 O , 6 c = ^-k\ (3.63) 
The anisotropic stress of the neutrinos can be determined by solving Eq. (3.13): 

°» = -mU*-4«)- lX - (3-64) 

Inserting eq. (3.64) into Eq. (2.37) the peculiar velocity of the neutrino fluid can be 
determined to be 

b , 3 ob 



Finally from Eq. (2.38), using Eq. (3.64), the difference between the two longitudinal 
fluctuations of the metric can be derived: 

ifo = (f>o(l + iRu) ~ + |ob. (3.66) 
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The momentum constraint derived in Eq. (3.4) is satisfied by virtue of the following identity 

a B k 2 7]-^ + n h — ^— = 0. (3.67) 



Notice also that the contribution of the Lorentz force is suppressed, for k — > with respect 
to Ob since, roughly, Fb ~ k 2 Q,B- 



4 Synchronous gauge description 

In the longitudinal description, the gauge freedom is completely specified and, hence, the 
so-called gauge modes are absent 4 . As a consequence [41], in the longitudinal gauge the 
metric fluctuations correspond directly to a particular set of gauge-invariant combinations, 
namely the Bardeen potentials [3]. Having said this, it is important to appreciate that 
longitudinal and synchronous descriptions should be regarded as complementary and not as 
opposite. The synchronous description, for instance, is more convenient when modes based 
on anisotropic stresses are discussed, as in the case of the present analysis. Furthermore, the 
known (old) problem of spurious gauge modes is completely settled by now since the work 
of Press and Vishniac [5]. In some specific calculations it turns out to be useful to profit 
from the extra gauge freedom inherent in the synchronous description . There are examples 
in the literature of isocurvature modes that are divergent, at early times, in the longitudinal 
gauge, but which are perfectly physical and regular in the synchronous coordinate system 
[37]. Finally, it is important to consider the formulation of the problem of initial conditions 
in the synchronous gauge, since various numerical codes solving the Boltzmann hierarchy 
use, indeed, the synchronous description. 

In the synchronous description the line element can be consistently perturbed in the 
form: 

ds 2 = a 2 (r])dr] 2 - a 2 (??)(<% - hij)dx % dx\ (4.1) 
where the perturbed element of the metric is given by 

S gij = a 2 hij, 5g 11 = — (4.2) 

Sometimes the parametrization Sg^ = 2a 2 (V> s <% + didjE s ) is also employed [41]. The 
connection between the two parametrizations is immediate by separating, in Fourier space, 
the trace of the perturbation from its traceless part: 



kikjh(k, rj) + 6£(/c, r?) ( kjtj - ^5ij 



(4.3) 



4 According to the usual nomenclature, gauge modes are those modes arising when gauge freedom is not 
completely fixed, as in the synchronous gauge. 
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where k l = k l /\k\. Clearly, by performing an infinitesimal gauge transformation the Sgoo 
and <5goi parts of the metric (which are not perturbed in the synchronous parametrization) 
also transform. The gauge modes can be exactly identified by requiring that Sgoo = and 
$9ij = 0) f° r g au g e transformations that preserve the synchronous nature of the coordinate 
system. The first gauge mode corresponds to a spatial reparametrization of the constant- 
time hypersurfaces. As a consequence, in this mode the metric perturbation is constant 
and the matter density unperturbed. The second gauge mode corresponds to a spatially 
dependent shift in the time direction. 

With the notation of Eqs. (4.2) and (4.3), the perturbed Einstein equations become 



k 2 i - ^h> = 

s 2 


-n 2 [Rj v + (i - R U )5 J + n B + n h s h + n c s c \, 


(4.4) 


k 2 ^ = -\n 2 


-7— — + UruOu + (i - Ru)0 1 ) + VbO h + n c e c 

4irap 6 


(4.5) 


h" + mti - ik 2 i = m 2 


l -(R u 5 u + (1 - R^) + ^ 


•> 


(4.6) 


h" + 6£" + 2Hti + 12H£' 


- 2k 2 i = m 2 \a B + \r v o v 




(4.7) 



where the background equations (2.3)-(2.5) have already been used to eliminate the energy 
and pressure densities 5 . By combining appropriately Eqs. (4.4) and (4.6), it is possible to 
obtain a further useful equation 

h" + Kb! = m 2 [2R y 5 v + 2(1 - R^ + 2Q B + n h 5 h + Q C 5 C }. (4.8) 

The evolution equations of the peculiar velocities and density contrasts of the various species 
of the plasma can be obtained by perturbing the covariant conservation of the energy- 
momentum tensor. The result is the following: 



s ' v = -\ e » + \ti> 


(4.9) 


4 2 

^ = --e 1 + -h', 


(4.10) 


*b = -*b + \ti, 


(4.11) 


^ = -^ + 7' 


(4.12) 


k 2 

6' u = -k 2 a v + —5 U , 


(4.13) 


k 2 


(4.14) 




(4.15) 



5 In this section the density contrasts and the peculiar velocity field are named in the same way as in the 
longitudinal gauge. It is understood that the density contrasts and the peculiar velocity fields are not equal 
in the two gauges and are related by the transformations listed below (see Eqs. (4.26) and (4.28)). 
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-W0c 



(4.16) 
(4.17) 



As anticipated the synchronous gauge modes can be made harmless by eliminating the 
constant solution for h and by fixing, for instance, the CDM velocity field to zero. These 
two requirements specify the coordinate system completely. 

The force-free solution can also be obtained within the synchronous description. In 
particular, for the metric fluctuations the solution is the following 



h 



To " G 



5 

k 2 rf. 



5 + 4fi„ 
6(15 + AR V 



-C 



60 



B 



, 2 2 

k t] , 



(4.18) 
(4.19) 



The constant introduced in Eqs. (4.18) and (4.19) has been defined in order to match the 
standard notation usually employed in the literature (see for instance [6]) to characterize 
the adiabatic (inflationary) mode. By solving Eqs. (4.9)-(4.17) the solution for the density 
contrasts 



5 C = 

and for the peculiar velocities 



5j = — $7r + 

8 V = —Ob + 
3 



2R u + 2, 



n B - -C 
30 3 



2R* + 3R v -5 n 2 
30i^ " 3 C 



, 2 2 

k f] , 



--^B + 

4 

- 7 n B + 

4 



R, 



l 



- — n B -- 

20 2 

^-ln B - c - 

20 2 



kW, 



e 7 = -^ v + 

4 



2i2„ + 3_ C 

sZb 

360 18 



2R 2 + 7R U — 
360^ 



-0 B 



23 + AR V 
15 + AR V 



(4.20) 
(4.21) 
(4.22) 
(4.23) 

(4.24) 
(4.25) 



can be obtained. The solution derived in the previous equations of the present section can 
be transformed into the longitudinal gauge by using the appropriate transformation, i. e. 



<pL = -^m + h)"+H(^ + hy], 



(4.26) 
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where the subscripts refer to the quantities evaluated either in the longitudinal or in the 
synchronous gauge. By setting = in Eqs. (4.18)-(4.19) and (4.20) -(4.25), we recover 
the standard adiabatic mode discussed in the appendix with the constant C determined as 

C = ^*>, (4.27) 

where <fio is the value of the adiabatic mode in the longitudinal gauge discussed previously. 
Of course, the solutions of the longitudinal gauge can also be directly transformed into the 
synchronous gauge using the appropriate transformation, which can be easily derived: 



£ = -*1>L - — / a(j) L dr], 
a J 

h = 6ip L + 6- J cuf> L dri, -2k 2 J J a(r/)<f> L (r/)dr]' 

6 (X) = g(\) + m ^ w + i)I /* ^ Lad ^ 

Qj J 

7 2 

9 (\) = 9 w _ !l_ r a(j)Ldr] _ (4 28) 

(X J 

Up to now we always studied adiabatic (or quasi-adiabatic) solutions. However, also 
isocurvature solutions can be generalized to include fully inhomogeneous magnetic fields. 
In order to look for isocurvature solutions it is useful to recall that in looking for perturbative 
solutions (both in the longitudinal and in the synchronous gauge) there are various small 
parameters. One is certainly kr), which is small outside the horizon. However, also 
and O c are small parameters deep within the radiation-dominated epoch. Let us make this 
statement more precise by considering the scale factor 

interpolating between the radiation-dominated phase for rj <C rji and the matter dominated 
epoch for rj S> In this case we can also write 

ttb,c = n b , c ^^-, (4.30) 
a(i]) + 1 

n B =Tl B — i— , (4.31) 
a(f]) + 1 

where the subscripts in Eq. (4.30) refer either to baryons or to CDM and where we took, 
for simplicity, rji = 1. 

We can then insert Eqs. (4.29) and (4.30)-(4.31) into the evolution equations for the 
perturbations in the synchronous gauge, regarding f2b,c as small paramters deep within the 
radiation-dominated epoch (rj — ► 0). The result of this procedure is summarized by the 
following solutions, which are valid in the case <tb = ^b/3: 

h ~ (-4^7/ + 6H b r/ 2 ), (4.32) 

Z~lli h r)-Ti h T?, (4.33) 
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for the metric perturbations and 



S 1 ~ y-^bv + ^bV 2 j , -n B (i -v + v 3 ), 

8 h = (1 - 2H b + 3n b r? 2 ), 

<5„ ~ ^--n b r? + 4n b r? 2 ^, 
5 C ~ 2n b + 3n b r? 2 , 



(4.34) 



(4.35) 
(4.36) 



(4.37) 



# 7 ~ -^ b fcV - ^ B (4r7 - 2r/ 2 + r/ 4 ), 
0„ ~ -^W - ^3(4^ - 2r/ 2 + ry 4 )|l 

#c = 0, 



(4.38) 



(4.39) 
(4.40) 




(4.41) 



for the fluid quantities. Clearly, for this mode, the adiabaticity condition is not satisfied. 
Furthermore, by transforming the solution to the Newtonian gauge, it is easy to check that 
the longitudinal fluctuations of the metric vanish for r\ — ► 0. In the limit J7b — ► this mode 
reduces to the baryon isocurvature mode already discussed in [37]. A similar solution can 
be obtained by changing f2 b ^c- 

5 Tight-coupling expansion in magnet oactive plasmas 

One of the analytical tools often employed in the theory of the CMB anisotropies is the 
so-called tight-coupling expansion [43]. Defining the differential optical depth as done in 
Eq. (2.28) the exact tight-coupling limit is realized when cjt — » oo and 1/r' — > 0. If tight 
coupling is exact, photons and baryons are synchronized so well that the photon phase-space 
distribution is isotropic in the baryon rest frame. Since the photon distribution is isotropic, 
the resulting radiation is not polarized. However, as the decoupling time approaches there 
is a regime where |l/r'| < 1 without being zero. The idea is then to tailor a systematic 
expansion in powers of |l/r'| and to consider not only the zeroth order but also higher 
orders depending upon the accuracy required by the problem. The evolution equations to be 
discussed are essentially the radiative transfer equations which are, in turn, the Boltzmann 
equations written in terms of the brightness perturbations. Within the notation employed 
in the present paper, these equations are derived in the equations (B.26)-(B.28) of the 
appendix. The brightness perturbations should be studied for each of the four Stokes 
parameters. The brightness perturbations Aj, connected with the first Stokes parameter I, 
are the temperature fluctuations, i.e. fluctuations in the intensity of the radiation field. As 
discussed in the appendix, recalling that n is the direction of the comoving-three momentum 



23 



of the photon and defining fi = k ■ n, the evolution of the brightness perturbation is 



-Ai + Aio + /it* - ^Ffe(M)<Sb 



(5.1) 



A[ + ikfiAi = (ip' - ik^(j)) + t 

where, in our notation, = 9^/{ik) and 

S* = Ai 2 + A Q0 + Aq2, (5.2) 
where, with obvious notation (discussed in detail in the appendix) 

Ai(fc, n, n) = £H) £ (2^ + !) a m(*. v)Pe<J*), 

i 

A Q (k,n, V ) = ^(-i) e (2e + l)A Q t(k,T,)Pe(»), (5.3) 

i 

An and Aq^ being the £-th multipole of the brightness function Ai and Aq . 

The function P2(a*) = (3/i 2 — l)/2 in Eq. (5.2) is the Legendre polynomial of second 
order, which appears in the collision operator of the Boltzmann equation for the photons 
due to the directional nature of Thompson scattering. The evolution equations for the 
brightness perturbations connected with the Q and U Stokes parameters are 

Aq + ikfiA Q = r'[-A Q + 1(1 - P 2 (ji))S ], (5.4) 
A'u + iknAu = t' - Au. (5.5) 

Finally, in order to close the system, the baryon evolution equation should be added using 
the definition = 0\,/{ik) either in Eq. (2.26) or in (3.8).: 

v' h + Hv h = -ik(J) - -[3iA n + v h ] - iW B (k), (5.6) 
a 

where 

w ^ = TO' (5 ' 7) 

a= Aj 1 - (5 - 8) 

In (5.6) the dipole of the brightness function appears directly. In fact, using the results 
derived in the appendix, the following chain of identities holds: 

7 = jj/c.F 7 i = 3An. (5.9) 

where jF 7l is the dipole of the photon phase-space distribution and An is the monopole of 
the brightness perturbation. 

The idea is now to expand Eqs. (5.1) and (5.4) in powers of the small parameter 
e = |1/t'|, i.e. the inverse of the differential optical depth of Thompson scattering. Before 
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doing the expansion, it is useful to derive the hierarchy for the brightness functions in full 
analogy with what is discussed in the appendix for the case of the neutrino phase-space 
distribution. To this aim, each side of Eqs. (5.1)-(5.4) and (5.6) will be multiplied by the 
various Legendre polynomials and the integration over [i will be performed. Noticing that, 
from the orthonormality relation for Legendre polynomials (Eq. (B.10) of the appendix), 



(5.10) 



and recalling that 



P (m) = 1, PM=^ W = - 1), P^) = \(^-^), (5.11) 

Eqs. (5.1)-(5.6) allow the determination of the first three sets of equations for the hierarchy 
of the brightness. More specifically, multiplying Eqs. (5.1)-(5.4) and (5.6) by Po(n) and 
integrating over /i, the following relations can be obtained 



A' w + fcAn = ip', 

A Q0 + AA Qi = — [Aq2 + A H - A Q0 ], 



(5.12) 
(5.13) 

(5.14) 



v ' + Hv h = -ikcp (3iAn + v b ) - -W B (k). 

a a 

If Eqs. (5.1)-(5.4) and (5.6) are multiplied by Pi(fi), the integration over /j, of the various 
terms implies 



-A' n - hA l2 + ^A M = - 1 y> + 7 > 
-A' Q1 - |/cA Q2 + |a qo = t'Aqi 



All + TTT-Ub 
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v' + Hv b = -ikcj) - — (3iAn + v b ) - -W B (k). 

a a 



The same procedure, using P2(m), leads to 

3 2 
-AJ 2 - -fcA B + -A; An = r' 



^A I2 -^(A Q o + A Q2 ) 



-A' Q2 - ^/cA Q3 + |fcA Q i = T' 



9 1 

Y^A Q2 - ^( A Q0 + Ai 2 ) 



v' b + Hv b = -ikcj) ( 3iA n + v b ) - -W B {k) 



a \ 



(5.15) 
(5.16) 
(5.17) 

(5.18) 
(5.19) 
(5.20) 



For I > 3 the hierarchy of the brightness can be determined in general terms by using the 
recurrence relation for the Legendre polynomials reported in Eq. (B.ll) of the appendix: 



k 



A « + t ' Ak = 2 7TT [Ml( *- 1) + 1)Al ^ +1 ) ] ' 

A' Qe + r'A Qe = ^[Mq^d -(£ + l)A Q(m) ]. 



(5.21) 
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We are now ready to compute the evolution of the various terms to a given order in the 
tight-coupling expansion parameter e = |1/V|. After expanding the various moments of the 
brightness function and the velocity field in e 

Ak = Am + eSu, 

AQ£ = Aq£ + €5q£, 

Vb = v h + e5v h , (5.22) 

the obtained expressions can be inserted into Eqs. (5.12)-(5.17) and the evolution of the 
various moments of the brightness function can be found order by order. 

To zeroth order in the tight-coupling approximation, the evolution equation for the 
baryon velocity field, i.e. Eq. (5.14), leads to: 

v b = -3iAn, (5.23) 

while Eqs. (5.13) and (5.16) lead, respectively, to 

Aqo = A K + Aq 2 . 

Aqi = 0. (5.24) 

Finally Eqs. (5.18) and (5.19) imply 

9Ai2 = Aqo + A Q2 , 

9A Q2 = A Q0 + A I2 . (5.25) 

Taking together the four conditions expressed by Eqs. (5.24) and (5.25) we have, to zeroth 
order in the tight-coupling approximation: 

A Q ^ = 0, I > 0, 

A u = 0, I > 2. (5.26) 

Hence, to zeroth order in the tight coupling, the relevant equations are 

v b = -Si A n , (5.27) 
A^o + fcAn = V' (5.28) 

This means, as anticipated, that to zeroth order in the tight-coupling expansion the CMB 
is not polarized since the quadrupole moment of the brightness is vanishing. Summing up 
Eqs. (5.15) and (5.17) and using Eq. (5.27) in order to eliminate from the obtained 
expression, we get to the following equation 

(a + l)A' n + HaA n - ^A I0 = ^(a + 1)0 + ^W B (k). (5.29) 
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Finally, the dipole term can be eliminated from Eq. (5.29) using Eq. (5.28). By doing so, 
Eq. (5.29) leads to the wanted decoupled equation for the monopole: 



-a-" a ' -ir' k 2 - r 

a + 1 o(a + 1) 



a + 1 3 



+ 1 12vrp b a 4 ' 



(5.30) 



In the limit Fb — > the usual decoupled equation for the evolution of the monopole is 
recovered. 

The presence of the magnetic field modifies the evolution of the monopole. Hence, also 
the dipole will be modified since the relation (5.28) stipulates that the monopole is a source 
of the dipole. 

To first order in the tight-coupling limit, the relevant equations can be obtained by 
keeping all terms of order e and by using the first-order relations to simplify the expressions. 
From Eq. (5.16) the condition Sqi = can be derived; from Eqs. (5.13) and (5.18)-(5.19), 
the following remaining conditions are obtained respectively: 

Sqo + <5i2 + 5 Q2 = 0, (5.31) 
^i2-^Q0 + 5Q2] = ^An, (5.32) 

^Q2--^Qo + <fo]=0. (5.33) 
Finally from Eqs. (5.31)-(5.33) the remaining relations are: 

£qo = |*I2, (5.34) 

5q2 = 4<5i2, (5.35) 

fa = ^A n . (5.36) 

Condition (5.36) can be also written 

Ah = ed l2 = 4^11- (5-37) 
15 T 

Now, from Eqs. (5.34) and (5.35), the quadrupole moment of Aq is proportional to the 
quadrupole of Ai, which is, in turn, proportional to the dipole evaluated to first order in 
e. But Aq measures exactly the degree of linear polarization of the radiation field. So, to 
first order in the tight-coupling expansion, the CMB is linearly polarized. Furthermore, as 
discussed in Eq. (5.30) the presence of the magnetic field modifies the evolution of both 
the monopole and the dipole. Thus, according to Eq. (5.37), also the polarization will be 
modified. 

All the the solutions discussed in the previous sections can be used to determine the 
initial conditions for the evolution of the brightness using the tight-coupling expansion. 
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Consider, as a possible example, the case when initial conditions are set deep within the 
radiation epoch. In this case, recalling Eq. (5.8), Eq. (5.30) can be written for, a <C 1, as 

A' ' + ^Ao = -co 2 ^ + V) - (5.38) 

107T/9 7 

where Ao = Aio — ip and cu = fc/\/3; we also posit a 4 p 7 = p^. The solution of Eq. (5.38) can 
be easily obtained in terms of arbitrary integration constants. These constants are fixed 
by specifying the initial conditions for the velocity field i>b, which determines, according to 
Eq. (5.23), the value of the dipole. Suppose, for instance, to be interested in the case of 
quasi-adiabatic initial conditions in the presence of the reduced Lorentz force. Then, from 
Eqs. (3.59) and (5.23), recalling that #b = ikvb, 

An = — ?? + — 3. (5.39) 

From Eq. (5.17) the initial condition for Ao can be determined to be: 

Aq^-^I + ^^o + ^^b-^b. (5.40) 

The same strategy can be applied to more specific cases, such as the one where the scale 
factor interpolates between a radiation-dominated phase and a matter-dominated phase, as 
discussed in Eq. (4.29). In this case the solution of Eq. (5.30) will be more complicated 
but always analytically tractable. 



6 Concluding remarks 

In this paper a systematic treatment of scalar perturbations has been discussed in the pres- 
ence of a fully inhomogeneous magnetic field. No specific configuration has been assumed 
and the results are, in this sense, rather general. Large-scale magnetic fields have been 
described in a fully consistent one-fluid MHD approach in curved space-time, which is par- 
ticularly suitable for the analysis of the low-frequency part of the plasma spectrum. In this 
approach the charged components of the plasma (baryons and electrons) are in thermal 
equilibrium at a common temperature. The neutral components of the plasma (photons, 
neutrinos and CDM particles) are not affected directly by the presence of the magnetic 
fields. However, since magnetic fields gravitate and appear in the perturbed Einstein equa- 
tions, also the initial conditions for the evolution of the neutral species are modified in a 
specific fashion. Since in MHD the Ohmic current is solenoidal, the baryon evolution equa- 
tion is affected by a reduced Lorentz force term whose characteristic form is well known in 
flat-space MHD. The combination of these different effects leads to a set of initial condi- 
tions for the CMB anisotropics, which is rather different from the standard adiabatic (or 
isocurvature) modes. 

The main new results of the present analysis can be summarized as follows: 
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• the problem of initial conditions for magnetized CMB anisotropies has been solved 
both in the conformally Newtonian gauge (more useful for theoretical calculations) 
and in the synchronous gauge (more appropriate for numerical discussions); 

• if the curvature fluctuations are adiabatic, magnetic fields modify the conventional 
adiabatic mode in a computable way and, as a consequence, the whole system of 
initial conditions for CMB anisotropies is modified; 

• for the modified adiabatic mode, deep outside the horizon, the density contrasts still 
satisfy the adiabaticity condition, however, as the horizon is crossed, a small non- 
adiabatic component develops; 

• if the fluctuations are isocurvature from the beginning, magnetic fields modify quan- 
titatively their amplitude, as shown in the case of the baryon isocurvature mode; 

• the tight-coupling expansion has been revisited in the presence of a fully inhomoge- 
neous magnetic field; 

• it has been shown that, because of the reduced Lorentz force, both the zeroth and 
first order in the tight-coupling expansion are modified and this allows the initial 
conditions for the evolution of the brightness functions to be computed reliably in 
both the adiabatic and isocurvature cases when fully inhomogeneous magnetic fields 
are present. 

It is appropriate to conclude with some remarks concerning the interplay between large-scale 
magnetic fields and CMB physics. In various investigations, limits on the magnetic field 
intensity are obtained on the basis of CMB considerations. Most of the times these limits 
refer to specific configurations. These discussions are certainly valuable; however, the limits 
obtained should always be associated with a specific set of analytical initial conditions. As 
far as scalar fluctuations are concerned, this was not done up to now. An example is Ref. 
[8] where no specific analysis of the initial condition problem was presented, but rather 
stringent limits were claimed. 

We must bear in mind that the analysis of CMB anisotropies requires the determination 
of a sizeable set of parameters. If the magnetic field were the only unknown, the task 
would be easier. However, this is unfortunately not the case. Hence, even within a specific 
magnetic field configuration, the limits should always refer to the specific set of initial 
conditions assumed in the analysis, as is normally done in the case when magnetic fields are 
absent. We do hope that future analyses along these directions may take the results of the 
present paper as a useful format for the systematic investigation of the effects of large-scale 
magnetic fields on the CMB anisotropies. 
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A Boltzmann hierarchy 



As discussed in the bulk of the paper neutrinos are collisionless particles. In the following 
we are going to discuss the collisionless Boltzmann equation following the notation discussed 
in the present paper. The phase-space distribution of massless particles (massless neutrinos 
or photons) is defined to be 

f{x\ P i ,ri)dx 1 dx 2 dx 3 dP 1 dP 2 dP 3 = dAf, (B.l) 

giving the number of particles in a differential volume of phase space. The function appear- 
ing in Eq. (B.l) is scalar under canonical transformations, Pi is the conjugate momentum. 
To derive and use the Boltzmann equation in curved space-time it is more convenient to 
work directly with the modulus q and the direction rii of the comoving three-momentum qi 
(where qi = qrii with riiri 1 = 1). Denoting the longitudinal degrees of freedom of the per- 
turbed geometry as in Eq. (2.6) and adopting the metric signature of Eq. (2.2), the relation 
between the comoving three-momentum and the components of the conjugate momentum 
are Pi = - ip) and P = q(l + 4>), 

The total variation of the distribution function can be written as 

Df df dx i df df dq df dn l fdf\ , n . 



Drj drj dr\ dx l dq drj drii drj \drj J 



coll 



where the collisional term has been kept for later convenience, but it is zero in the case of 
neutrinos. The collisionless part of Eq. (B.2) can be perturbed around a configuration of 
thermal equilibrium by splitting the phase-space distribution as 

f(x\ Pj,ri) = f Q (q)[l + fM{x\q, n jt r?)], (B.3) 

where fo(q) is the unperturbed phase-space distribution, which only depends upon the 
comoving three-momentum. 

Inserting Eq. (B.3) into Eq. (B.2), only the terms that are first order will be kept. 
Since the unperturbed phase-space distribution only depends upon the comoving three- 
momentum, Eq. (B.2) becomes: 

^W^ + ^-n^-i^) , (B.4) 
drj dx l ding fo\dr) ) co n 

where the geodesic equation has been used in order to obtain the expression of the time 
derivative of q to first-order in the amplitude of the metric perturbations, i.e. 

= qi f/ - qm&f. (B.5) 

drj 

Up to this point the derivation is valid for both photons and massless neutrinos. 



30 



A.l Massless neutrinos and photons 

Going now to Fourier space and denning the reduced phase-space distribution for massless 
neutrinos as 

F v (k,n,ri) = — , B.6 

J q s dqf 

Eq. (B.4) becomes, in the absence of collision term, 



+ ik[iT v = 4(tp' - ikficj)), (B.7) 



where n = n-k. The factor 4 appearing in Eq. (B.7) follows from the explicit expression of 
the Fermi-Dirac phase-space distribution and observing 

The reduced phase-space distribution can be expanded in series of Legendre polynomials as 

T v % n, 7]) = ^2(-i) e (2£ + l)F vt (k, r})P t {n). (B.9) 

i 

Equation (B.9) will now be inserted into Eq. (B.7). The orthonormality relation for Leg- 
endre polynomials, 

£ Pe(»)PA») = (B.10) 
together with the well-known recurrence relation 

(* + l)P m 00 = (2£+l)tiPe(ti)-ePe-i(») (B.ll) 

allows to get a hierarchy of equations to be obtained for the various multipole moments. 
The procedure is to take the various moments of both sides of Eq. (B.7). In doing so, 
expressions like 



(B.12) 



will appear; they can be evaluated by using Eqs. (B.10) and (B.ll). The full form of the 
Boltzmann hierarchy is then 

^o = -fcF„i + V, (B.13) 
Ki = |[Jvo-2^v2] + |fe^ 

Ke = ^C- 1 ) + i^+i)]- ( B - 14 ) 
Equation (B.14) holds for I > 2. 
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The components of the energy-momentum tensor can be connected with the monopole 
and dipole of the distribution function 

d 3 P P„P 



^ p fV'Pj'V)- (B.15) 



Recalling the connection between conjugate momenta and comoving three-momenta we get, 
to zeroth order: 

Pu = -I d 3 qqf (q), (B.16) 

CI J 

while to first order the density contrast (monopole), the peculiar velocity field (dipole) and 
the quadrupole become 

1 



dUJ r l/ (k,h, 77) = Fvo, (B.17) 

Q ? ' r _ _ 3 

e » = T7T dto(k-n)F v (k,n,ri) = -KF vl , (B.18) 

l07T J 4 



{k-hf- 1 - 



F v {k,n,r,) = ^. (B.19) 



Inserting Eqs. (B.17) and (B.19) into Eqs. (B.13)-(B.14), the system following from the 
perturbation of the covariant conservation equations can be partially recovered 

%, = -p v + 4i/>', (B.20) 
k 2 

K = -jtv - k 2 a„ + k 2 <P, (B.21) 

< = ^ - AfcF^ (B.22) 

with the important addition of the quadrupole (appearing in Eq. (B.20)) and of the whole 
Eq. (B.22), which couples the quadrupole, the peculiar velocity field, and the octupole J- v %. 
Equation (B.22) is important. After neutrino decoupling, when initial conditions are set, 
Fv3 = 0. 

Equation (B.4) can also be made explicit in the case of photons. In order to do so the 
collision term should be specified. Before writing down the Boltzmann equation for the 
photons it is useful, in order to match with the standard notations, to define the photon 
brightness perturbation, conventionally denoted by A: 

f(x i ,q,n j , V ) = f (j^j, (B.23) 

where fo(q), i.e. the unperturbed phase-space distribution, denotes now the Bose-Einstein 
distribution. Comparing Eq. (B.3) with the expression given in Eq. (B.23) and expanding 
for A < 1, the two definitions are connected as 

} UlnJ ' Jq^dqfo ^ 
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where T 1 is the reduced phase-space distribution for the photons defined in full analogy with 
the case of the massless neutrinos, i.e. Eq. (B.6). The second identity of Eq. (B.24) follows 
from the first and from Eq. (B.8). Using the relation between the brightness perturbation 
and the phase-space distribution, the collisionless part of the Boltzmann equation can be 
written as 

A' + ikfi(A + (/))= ip'. (B.25) 

The brightness perturbation, as defined in Eq. (B.23), corresponds physically to a pertur- 
bation in the first Stokes parameter (conventionally denoted by I), i.e. to the sum of the 
squared amplitudes of the radiation field. By specifying the collision term, the full form of 
the Boltzmann equation for photons can be written as 



Aj + ikfiAi = {ip' - ikn4>) + r' 
A' Q + ikfiA Q = t' 



-Ai + A 10 + nv h - -P 2 (n)S 



(B.26) 



A Q + i(l-JM M ))So 



(B.27) 

A(j + ikfiA v = -t'Atj. (B.28) 



Concerning Eqs. (B.26)-(B.28) a few comments are in order. The brightness functions Aq 
and Au correspond, respectively, to the Stokes parameters Q and U. The parameter Q 
being the difference of the squares of the amplitudes of the radiation field, is sensitive to 
the linearly polarized radiation. Since Q and U (unlike I and V) change under rotations, 
once Q is included also U must follow. The source term appearing in Eq. (B.26), So is, 
within our conventions: 

6*0 = A K + A Q0 + Aq2, (B.29) 

where, with obvious notations, A# denotes the £-th multipole of Ai and similarly for Aq. 
The peculiar velocity field for the baryons has been written as 

v h = |. (B.30) 

The notation v\, is preferred here in view of the application to the tight-coupling expansion of 
the Boltzmann hierarchy for the photon brightness. Finally, r' = x e n e o^-^- is, as previously 
introduced, the differential optical depth for Thompson scattering. The evolution equation 
for «b can be easily obtained, for instance by Fourier transforming Eq. (4.26) and by using 
Eq. (B.30): 

v' h + Hv h + ik(f } + - f 3iAn + v b ) = - \ FB } k \ (B.31) 
a V / 47tk a 4 /9b 

having defined 

a = ^. (B.32) 

4 pr 

Equations (B.26)-(B.28) and (B.31) are discussed in Section 4. 
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